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A STEADY STATE LONGITUDINAL MANPOWER PLANNING MODEL 



WITH SEVERAL CLASSES BY MANPOWER 

by 

Richard C. Grinold 
ABSTRACT 

In the Navy's officer system there are several classes of 
manpower (e.g. pilots and submariners) that perform specialized 
jobs and also can perform several types of non-speciali zed jobs 
(e.g. military training, personnel management, etc.). This report 
is concerned with the general problem of allocating the different 
types of jobs among the several classes of manpower. 

The report describes a model that constructs a personnel 
inventory by rank for each of several manpower classes (pilots, 
etc.) and then allocates those people to the specialized and 
common jobs that they are allowed to do. 

The idea of the model is to allow a policy maker to quickly 
reconcile billet requirements with the reality of available acces- 
sions, job sharing targets between classes, and continuation rates 
of the different manpower classes. 

Different allocations can be produced either by assuming 
the values of a few key variables, or they can be generated using 
an optimization scheme that sets "allowable" percentage errors. 

Four optimization variants based on this idea are described. The 
report contains some typical data, the results of calculations, 
and a description of computer programs used to solve the problem. 



A STEADY STATE LONGITUDINAL MANPOWER PLANNING MODEL 



WITH SEVERAL CLASSES OF MANPOWER 



by 



Richard C. Grinold 



1 . Introduction 

This paper outlines the construction of a manpower plan- 
ning model for a system in which each of several classes of man- 
power are assigned to several categories of jobs. The jobs are 
either specialized in that only one type manpower can perform 
the job or general in that more than one manpower type is quali- 
fied to do the job. 

Both jobs and manpower are broken down by experience 
level. A job with a certain experience level must be filled with 
a person that has the same experience level. 

The model assumes a steady state. In each year the job 
requirements are the same. In each year the same number of people 
enter the system (bottom level entry only) , the same number leave 
the system from each experience level, and the same number move 
up from each one level so that the inventory of people at each 
exerpience level and in each manpower class remains the same . 

There are several items represented in the model: 
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(i) The requirements for jobs by experience levels. 

(ii) The requirements for people by experience levels. 

(iii) The job sharing targets for common jobs; the fraction 
of common jobs we would like to see allocated to each 
manpower class. 

(iv) The people sharing targets for manpower classes; the 
fraction of each manpower class that we would like to 
see allocated to each type of job. 

(v) The rate of accessions in each manpower class. 

(vi) The retention of manpower in each class; i .e . the 

relationship between accessions in each class and the 
inventory of people by experience level in each class. 

The model is designed to study the interaction of these 
factors. It is easy to see from the list above that items (i) 
through (iv) are interrelated and unless care is taken in specify- 
ing these goals they will be inconsistent. Also items (ii) , (v) 

and (vi) are strongly related in that (v) and (vi) determine the 
inventory of people in each class at each experience level. Only 
in the most fortunate circumstances will this be consistent with 
the personnel requirements (ii) . 

Items (i)-(iv) are targets. They may be unrealistic and 
inconsistent but they can be set by the manpower planner. Items 
(v) and (vi) , however, depend to some extent on the behavior of 
personnel in the system (vi) and on the system's ability to 
attract qualified people to each class (v) . The reader 
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should keep in mind this caveat and be prepared for shifts in 
personnel behavior in response to changes in other system parameters 
or in external factors . 

The model is a laboratory for testing the relationship 
between these factors. We can try one set of policies and look 
at the discrepancy between job requirements and actual require- 
ments as well as the discrepancy between personnel requirement 
and the actual inventory in each class. 

The model may be used in several ways. The easiest is for 
the planner to stipulate a policy and then examine its effects. 

A more difficult procedure is for the planner to stipulate a 
range of policies and then use some type of optimization scheme 
to select a policy within the range. This paper contains four 
examples of this type of optimization. Each is based on the notion 
of a penalty function the measures the discrepancy between desired 
job and personnel requirements and what is actually provided. 

We take two types of penalty functions : piecewise linear and 
quadratic, and we examine two broad sets of policies: allocating 
people to jobs using the job sharing rules ((iii) above), or 
allocating people to jobs using the people sharing rates ((iv) 
above) . This gives us four combinations and thus four distinct 
optimization models. The piecewise linear penalty functions lead 
to linear programming models and the quadratic penalty functions 
to the minimization of a quadratic form subject to linear equality 
constraints . 
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The model is motivated by a study of the U. S. Navy's 
officer corps. We shall carry an example using that system 
throughout the text in order to illustrate each idea. 

The model is based on a longitudinal manpower flow model. 

This type of system is described in depth in Grinold and Marshall 
[1 ] , Chapter 3. The use of piecewise linear penalty functions 
to measure the discrepancy between actual performance and stated 
objectives is commonly called "goal programming." This idea 
has been extensively developed by Charnes and Cooper; a good 
review can be found in [3], The use of quadratic penalty functions 
is quite common in the optimal control literature; its use in 
more behavioral settings was pioneered by Holt, Modigliani, Muth 
and Simon in [2]. 

The paper consists of several short sections each dedicated 
to a specific point and most illustrated by the example carried 
thoughout the text. The model's structure is given in Sections 1-14. 
Section 15 is a review that gathers all the definitions presented 
to that point. Sections 16 to 20 describe alternative ways to 
choose an allocation, the idea behind our use of penalty functions, 
and the specific construction of the piecewise linear and quad- 
ratic penalty functions. 

Appendix A describes the two linear programming models 
that arise from the use of piecewise linear penalties and 
Appendix B the two models that stem from the use of quadratic 
penalties. Appendix C contains some sample solutions. 
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The organization reflects the relative importance of 
the topics. The structure of the model is the most important 
and it is stressed and reviewed in Sections 2-14. The use of 
optimization to select an allocation is, in the main, merely a 
device to circumvent the difficulty of having a wide range of 
policy choice. Once the optimization rule is set the planner 
has a direct route from policy to result. It is the variation 
in results that comes from changes in policy that will be of 
most interest to the planner. The optimization is a device to 
help make that connection. 

One final point should be made before describing the 
model's structure. This is an aggregate planning model. It is 
intended to test policies that will, in turn, provide a foundation 
for the day-to-day operating of the system. This model certainly 
will not show us how to operate the system. No model can answer 
all questions simultaneously and ours is no exception. 
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2 . Manpower 

There are K different classes of manpower indexed 
by k = 1,2,..., K. The classification scheme is, of course, 
directly related to the objectives of the model builder. In 
general, the classification should be fine enough to capture 
the important substitution possibilities and economical in avoid 
ing the listing of all possibilities in a futile attempt to repl 
cate reality. In our example there are five manpower types. 



Manpower 

Type 


Officer 

Designator 


Description 


1 


llOx 


Women 


2 


lllx 

(and 116x) 


Surface warfare 


3 


112x 

(and 117x) 


Submarine warfare 


4 


1 31x 
and 
1 39x 


Pilots 


5 


132x 

and 

137x 


Naval flight officers 



Table 2-1: The manpower classes. 
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3 . Stages 



Each officer's career is broken down into I-stages indexed 
by i = 1 , 2 ,..., I. The stages can be defined in many ways. The 
simplest is by period of service; if we track officers for 26 
years then the index i would go from one to 26. In our example 
the stages roughly coincide with the time period in which officers 
hold a certain rank. In general for i = 1,2,..., I, stage i 
will run from time of service s^_^ to s . . The maximum length 
of service is Sj = S and s^ = 0 . 

A person in manpower class k and in stage i of their 
career will be called a type (i,k) person. 



Stage 


LOS 


Description 


1 


0-2 


Ensign 


2 


2-4 


Lieu tenant-- J . G. 


3 


4-9 


Lieutenant 


4 


9-14 


Lieutenant Commander 


5 


14-19 


Commander 


6 


19-26 


Captain 



Table 3-1. Stages or experience levels. 
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4 . Manpower Flows 



We assume a steady state model. In each year the same 
number of people will enter each manpower class, and the inven- 
tory of people in each year of service and in each manpower class 
will remain constant. 

For class k = 1,2,...,K let y^ be the number of people 
entering manpower class k, and let a(k,s) be the fraction of 
these accessions that remain in the system for s years. The 
index s runs from 0 to S years and the function a(k,s) is 
decreasing in s. The function a(k,s) is frequently called the 
survivor curve. 

Stage i of a person's career runs from length of service 

s i-l to s i (s 0 = * T ^ e nuitl ^ er °f people in stage k is 

w .,y, where w., is defined by 
ik-^k ik 



w 



ik 



= / a (k, s) ds 



*i-l 



In our example, let's take a hypothetical survivor curve 
for pilots (class 4) . We follow these officers for 26 years. 
Figure 4-1 shows the survivor curve. The shaded areas are 
integrals under this curve for the duration of each stage. 

The calculation of w^ k can be approximated using a 
discrete form of the survivor curve. 

The coefficients w^ can b e interpreted as the amount 
of time a person in manpower class k is expected to spend in 
stage i. A crude way to view this is to say that w^ is the 
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product of the probability of reaching stage i, and the length 



of stage i. 



The probability will be 



w . , / ( s . 
ik' 1 



s 



i-1 



) . 





i 


w . , 
lk 


s . - s . . 

l l-l 


Prob . of 
reaching 
stage i 


s 


i 


1.97 


2 


0 .985 


T 


2 


1.85 


2 


0.925 


A 


3 


2.50 


5 


0.5 


G 


4 


1.40 


5 


0.28 


E 


5 


1.22 


5 


0.25 


S 


6 


0.54 


7 


0 . 77 



Table 4-1. The expected waiting time. 
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5. Jobs 



There are J different types of jobs indexed by 
j = 1,2, ,J. In our example we consider seven job types: 



Job No. 


Billet 

Code 


Description 


1 


1000 


General, nonwarfare billets 


2 


10 50 


General, warfare billets 


3 


1110 

and 

1160 


Surface 


4 


1120 

and 

1170 


Subsurface billets 


5 


1310 

and 

1390 


Pilots 


6 


1320 

and 

1370 


Naval Flying Officer 


7 


1300 


General Aviation 



Table 5-1. Definition of Jobs 
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V£> 

CO 




V£> 

CM 



nr 

CM 



CM 

CM 



O 

CM 



LD 

CO 



*3* nr 
iH CO 



m 

co 



nr 



CM 

co 



co 



o 
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Figure 4-1. Length of service 



6 . Billet Requirements 

The I by J matrix B contains the billet requirements. 
Since the career stage is a measure of experience (e.g. length of 
service or rank), the requirements are broken down by stage. 




The billet requirements matrix for our example is shown below. 

The billets will be identified by the index pair (i,j). 
Thus we shall speak of type (i,j) billets and type (i,k) people. 



JOBS 







1 


2 


3 


4 


5 


6 


7 


s 


1 


409 


0 


2095 


688 


961 


521 


0 


T 


2 


1008 


0 


1883 


734 


1998 


1212 


1 


A 


3 


1806 


378 


2080 


84 4 


3572 


1399 


636 


G 


4 


1495 


599 


1464 


810 


1780 


547 


726 


E 


5 


1031 


470 


926 


501 


678 


55 


902 




6 


592 


490 


321 


149 


0 


0 


387 



Table 6. The billet requirements matrix B. 
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7. Billet Sharing Array 



For each i the K by J matrix F 1 gives the share 

(fraction) of the type (i,j) billets that should be performed 

by manpower type k. The elements of F 1 are denoted f^ . , 

K 3 

they are nonnegative, and the rows sum to one, ^ f^j = 1. 

If fj^j > 0 , we interpret the fraction as a goal; we would 
like that fraction of the (i,j) billets filled from manpower 
class k. However, we interpret fjV = 0 to mean that manpower 
class k is not qualified to fill billet (i,j). 

There is also possible confusion between the idea of job 
sharing and a people sharing concept. A manpower class may per- 
form two types of jobs: the jobs for which it is uniquely 
qualified and other jobs that are shared among several manpower 
classes. Another way to look at the allocation of the common 
jobs is to stipulate a fraction of the inventory of (i,k) 
people that should be assigned to job j. That idea will be 
considered in Section 12. 

To save space, we shall only give one of the matrices 
F 1 (for i = 3) used in our example. 



JOBS 




Table 7-1. The job sharing matrix F 1 for i = 3. 
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8 . Target Allocations 



The number of (i,k) people desired in job (i,j) is 



given by 



t J . = f l . b . . 

k] kj ij 



Recall that b^ is the (i,j) billet requirement and fjV is 
the fraction of those billets to be filled by manpower class k. 

The target allocations used in our example (for i =» 3) 
are shown in Table 8-1. 



JOBS 




Table 8-1. 



The target allocation 
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9 . Actual Allocation 

Let the variable a£ . give the actual number of (i,k) 

people assigned to billet (i,j). If it was possible we would 

like a^j = t^ ^ . That would ensure that all targets are met. 

That, however, is usually impossible. Indeed, reconciling 

the a^ and t^ is the purpose of this model. 

The variables a, 1 . do not mean that a. 1 . individuals 

kj k j 

from class k are locked into type (i,j) billets when they 
are in stage i. It does mean that at any time a^j people 
from class k are filling type (i,j) billets. Any particular 
(i,k) person, may spend stage i in several billets. 



If fj~ . = 0, then an (i,k) person is not qualified to 
kj 

fill an (i,j) billet. So if f£ . = 0 , we shall require that 

a. 1 . = 0 . 
kD 
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10 . 



Conservation 



Recall that is number of type (i,k) people. 

The allocation must satisfy the following conservation 
relation 

J i 

\ a kj = W ik^k for a11 (i ' k) 

1=1 

This says that all. manpower fills some job. 
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11. Excess and Deficit 



One of our main objectives in constructing this model 
is to compare the allocation and target. We could look at 
all the discrepancies a^ k - t^. That would be very general, 
but it may not be comprehensible. We have chosen to concen- 
trate on more aggregate measures of excess and deficit: the 
actual number of (i,k) people and the actual number of 
(i,j) billets filled as compared with the targets. Our 
notation for these measures is 



K 

x. . = y a. . the number of poeple assigned 

ID k =l k D ^ ^ 

to billet (i,j) 



x^j - b_^j the discrepancy between billet 

(i,j) assignment and requirements 



z 



ik 



P 



ik 




J 



l 



j = l 




the number of (i,k) people 



the total requirement for type 
(i,k) people 



Z ik ' Pik 



the discrepancy between the inventory 
of (i,k) people and requirements. 
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12 . People Sharing 



Recall that we defined our targets from the viewpoint 
of billet sharing. Given any allocation we can calculate the 
people sharing fractions. Let g 1 be the fraction of type 

kj 

(i,k) people that are assigned to job j, thus billet (i,j). 

The g^ . are given implicitly by the equation 

a ij = 9kj z ik 

As mentioned above, an alternative model could be constructed 
using the people sharing fractions gjV as input targets 
and the billet sharing fractions f^ as outputs. In that 
case we would start with people requirement and the 

people sharing rules g^ and define target allocations by 

b kj - g kj P ik" 

Then we can calculate b. . from 

ID 

b ij ” | g kj p ik* 

In any case, for consistency, the following relations should 

hold 
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(i) «kjpik = f ij b ij 



(ii) I gjt = 
j-1 k:l 



= 1 ==> p 



ik 



J 

1 f£ -b. . 
j=l k => ^ 



(iii) Ji = 1 ==> = Ji 9 



To assure consistency we should either start with f 
and b and then calculate g and p, or start with g and 
p and calculate f and b. If we have a desire to start with 
g and b, then one suggestion is to guess p, see if it works, 
and then do some revision. 

From the values of b and f (for all i) used in this 
paper, we can compute the personnel requirements, shown below. 



Manpower Class 

1 2 3 4 5 



s 


1 


286 


2156 


713 


981 


537 


T 


2 


70 6 


2034 


794 


2049 


1253 


A 


3 


1264 


2540 


1028 


4093 


1790 


G 


4 


1046 


1988 


1019 


2376 


991 


E 


5 


722 


1316 


657 


1328 


540 




6 


414 


655 


283 


338 


249 



Table 12-1. The personnel inventories, p^ 
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13. Stage Substitution 



The astute reader will notice that we have not allowed 
substitution between stages. Thus a type (i,k) person is not 
allowed to fill a type (i+1, j) billet or a type (i-1, j) 
billet. 

Our formuation actually allows such assignments, however, 
we are assuming that they net out to zero. That is for each 
(i,k) person filling a stage i+1 billet, there is an (i+l,k) 
person filling a stage i billet. 

One way to study the possibility of stage substitution 
is to vary the length of the intervals (s^_^, s^) that 
define the stages. However, the billet requirement are pre- 
sumably set with some experience level in mind. Therefore 
major shifts in the definitions of stages would require adjust- 
ments in the billet requirements. 
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14. Cost 



The cost of any allocation can be calculated using 



I l c ik w ik y k • 



where c^ is the annual cost of an (i,k) person and = z 

is the number of (i,k) people. The cost data c^ should 
include salary, benefits, retirement; training, promotion, and 
recruitment costs. Notice we can define c^ as ^ c^ ^ , 
and then rewrite the cost as ^ c k^k' We interpret c^ as 
the average career cost of each accession into manpower 
class k. 
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15. Review 



We have constructed our model. The remainder of the 
paper is devoted to ways of calculating particular assign- 
ments and contains a sample calculation. 

We summarize here the notations and definitions pre- 
sented to this point. 



manpower class 


indexed by k = 1,2,...,K 


stages of a career 


indexed by i = 1,2,..., I 
defined by length of service 
(s^, s ^ / ... , s ^ ) 


jobs 


indexed by j = 1,2, . . . , J 


people 


indexed by stage and class (i,k) 


accessions 


y^, the number of people entering 
class k per year 


stage inventory 


w ik^k -*- s t ^ ie num ^ er °f (i,k) 
people; y^ is the accession rate 

and w., is obtained from the 

ik 

survivor curve and the stage 
definition. Note z^ = w ij c ^] c ' 


billets requirements 


indexed by stage and job (i,j); 
b^j is the billet requirements 


billet sharing 


indexed by (i,j,k); f£j is the 
fraction of (i,j) billets that 
should be filled from class k. 
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target allocation 


tjV . The number of (i,k) people 
desired in billet (i,j) 


actual allocation 


a^ . the number of (i,k) people 
K D 

in billet (i,j) 


class sums 


x. . and b. . are respectively 
the number assigned to billet (i,j) 
and the requirement for billet 

(i, j ) 


job sums 


z.^ and Pjj^. are respectively 
the number of (i,k) people, 
and the target for (i,k) people. 


people sharing 


g^ j , the fraction of (i,k) people 
assigned to billet (i,j); 

a kj = Z ik' 


cost 


c., is the annual cost of an 

xk 

(i,k) person. The cost of a particular 
allocation is ^ c ik z ik* 
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16 . Choice of an Allocation 

Up to this point we have set out a model structure. 

The next task is to devise one or two more procedures for 
selecting allocations. There are two general approaches to 
this task. 

On approach is to use an ad hoc rule. For example, 

we could fix accessions at projected rates for each class. 

That would fix the variables and thus z ^ = w^y^- 

Next we could specify the people sharing rules g* . . That 

K 1 

would give us the allocation 

The second approach is to use some type of optimization 
to select a policy. This optimization is either based on a 
trade-off between the cost of any allocation and some measure 
of its quality, or simply a measure of the quality of the 
allocation. 
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1 7 . Quality of an Allocation 



In order to choose an allocation using optimization 

we need some measure of the quality of that allocation. 

In our model we have already decided to focus on the 

discrepancy between people inventory and targets (z^ k ~ 

and the discrepancy between billet assignments and requirements 

(x. . - b . . ) . 
il ID 

We have selected the simplest form of quality measure 
that makes sense. It is a penalty function that is zero for 
a perfect allocation and positive for others. Thus it measures 
the lack of quality. The penalty function has three properties. 

(1) The penalty can be written as 



l 




l . . (x . . ) 

ID ID 



where 

(2) h ik (z ik ) and ^(x^) are nonnegative, convex and equal 
to zero if z ±k = p ik and x^ = b^ respectively. 

(3) Parameters 6^ and 4^ measure a unit upper and lower 
percentage error in meeting billet requirements 






Similar parameters 
cies between actual 
ments . 



if 


x . . 


* 

•I 

rQ 

CD 

+ 
i — 1 

I! 




ID 


ID ID 


or 


x. . 


= (1 - \lf. -)b.. . 




ID 


ID 1 D 


ik 


and 


6 .. are used f 

ik 



Inventories of people and require- 
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h iK U ik> = 1 



if z. k - <1 + 4> ik >P lk - 



°r z ik - (1 - « ik )P ik • 



Item (1) above says that the quality measure is 
separable, item ( 2 ) that is nonnegative, convex, and zero if 
the assignment is exactly on target. The third item requires 
more discussion. Item (3) is designed to answer the question: 
How do we compare a 4 percent shortfall in meeting a critical 
target with a 10 percent short fall in meeting a less critical 
target? Our answer to this question is to take a single 
target as a benchmark and to say arbitrarily for that target 
that a certain percent over and a certain percent under the 
target yields an error of one. Then for any other target 
we can compare the seriousness of deviations with our benchmark. 
We say a fraction 0^ over target or under target is 

as serious as the deviations we have established for our 
benchmark . 

For example, we could take stage 3 pilots as our bench- 
mark (i = 3, j = 5) , and take the percentage under as 4% 

( 4*25 = 0*04) and percentage over as 10 % ^35 = 0 . 10 ) as the 
unit serious deviation for our benchmark. Now take any other 
category, for example, stage 4 general warfare billets 
(i = 4, j = 2) . Then we ask how much under target would the 
assignment to these billets have to be in order to be as 
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serious as a 4% shortage of stage 3 pilots. In this way we 



can try to make the essential judgments about the 
between different categories. 

We shall now give two examples of penalty 
and then give some practical examples of how such 
could be used. 



trade-offs 

functions 
a criterion 
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18. Piecewise Convex Penalties 



One of the simplest ways to construct the penalty measure 
for deviations from target is to use a piecewise convex function. 
To simplify notation we should consider the case of a discrepancy 
between the inventory of people and desired inventories and 
we shall drop the (i,k) subscripts. 

We need a piecewise linear convex function h(z) that 
satisfies 



(i) 


h (z) = 0 


if 


x = P 


(ii) 


h(z) = 1 


if 


x = (1 + <J>) p 


( iii) 


h (z) =1 


if 


x = (1 - 6 )p 



Such a function is shown in Figure 2. This function 
can be represented in several ways. One of the simplest is as 
the maximum of two linear functions 



h (z) 



max 



1 (ZIE) 
4) 1 p ’ 




n 



This functional form is flexible and easy to manipulate. 
We shall see below, that it can be used to obtain an allocation 
using a linear program. Two examples of this are contained 
in the appendix. 
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Figure 18-1. 



Piecewise linear convex function. 
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19 . Quadratic Penalty Functions 

Quadratic penalty functions are also quite easy to use. 
Unfortunately, we cannot get a quadratic function to satisfy 
all of our conditions. For example let us write 



We want to have 

(i) £ (x) = 0 if x = b 

(ii) £ (x) = 1 if x = (1 + 6)b 

(iii) £ (x) =1 if x = (1 - ^)b 

(iv) £ (x) 0 

where item (iv) implies £ (x) has its minimum at x = b. 

These four conditions cannot be met by a quadratic 
function unless Q = ip ; that is unless there is a symmetry 
between being under and over target. In the symmetric case 
we can write £(x) as 



When the penalties are not symmetric, we must relax one of 
our four conditions. If we relax condition (iv) then £(x) 
becomes 



£(x) = a(x - b)^ + 3 (x - b) + y. 




when 0 = if>. 
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h (x) 



1 ( x - b\ 2 1 (e - ip) / x - b x 

6ip \ b / “ 2 dip { b ’ * 

This function has its minimum midway between (1 + 0)b and 
(1 - ^)b and it has a value of 

(0 - 4>) 2 

2dip 

at the minimum. 

This approximation is useful if 9 and ip are similar. 
However, the approximation becomes much worse when 9 and ^ 
are much different. 

Two quadratic penalty models are described in the 
appendix . 
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20 . Cost Quality Trade-Off 



We saw above that any particular allocation has an 
annual cost and also a penalty associated with its deviation 
from targets for billet requirements and personnel inventories. 
By placing different weights on the cost and penalty we can 
obtain a family of objectives that will lead to efficient 
allocations . 

The cost-penalty frontier is shown below. Notice 
the minimum cost solution is zero, since we would not have 
any accessions; i.e. y = 0, thus no cost. 

We write our objective as 




As the parameter X increases from 0 to 1 more emphasis is 
placed on the cost and less on the penalty. At X = 0 , we 
minimize the penalty and ignore the cost; at X = 1 we 
minimize cost and ignore the penalty. For intermediate values 
of X we establish a trade-off between cost and penalty. 

This curve should be used with some caution since the costs 
are expressed in real units (dollars per year) , but the 
penalty is expressed in ad-hoc units. 
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penalty 



minimum 

penalty 



APPENDIX A 



LINEAR PROGRAMMING MODELS 



This appendix describes two linear programming approaches 
for selecting an allocation. The discussion is brief and is 
intended for those familiar with linear programming. The first 
model is called LPB for billet sharing and the second is called 
LPP for people sharing. 

In both models the idea is to choose the number of 
people assigned to billet (i,j), and the number of accessions 
in manpower class k; these are respectively denoted and 

y^. The allocation a^V is then fixed using an ad hoc rule. 

We have selected two rules and thus get two different linear 
programming models . 

The first rule is for sharing billets 
Rule B (billet share) 



(A-l) 





x 



ij 



Recall that f ^ . is the fraction of type (i,j) billets that 
we desire to have satisfied from manpower class k. Rule B 
allocates the error in meeting billet requirements among man- 
power classes so that there is a constant percentage error in 
each class's allocation. If there is a percentage error of 
H • • in meeting the billet (i,j) requirement, i.e. 
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then the billet sharing rule means that the percentage error 

in the assignment of type k manpower to billet (i,j) will 

also be n . . . 

ID 

The people sharing rule is 



Rule P 



(A-2 ) 



a kj ^kj Z ik 



where z ^ is the number of type (i,k) 
type (i,k) people are allocated among different 
people sharing parameters gjV . if there is an 
in the inventory of type (i,k) 



people. The 

jobs using the 

error of n . . 

lk 



'ik 



(z ik ~ Pik 1 



then the percentage error in the assignment of type (i,k) 
people to job j will be r)^ (recall tjV = g^p^) • The 
billet sharing rule leads to a larger linear program than the 
people sharing rule. 

The constraints of the linear program for billet 
sharing are given below. The objective is the same for both 
programs; it will be discussed after we examine the constraints. 
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Table A-l. Constraints for LP- 



The first set of constraints is simply the conservation 

of people. There will be w^y^ type (i/k) people and under 

the billet sharing rule f^.x. • of them will be assigned to 

X 1 11 

billet (i , j ) . 

The second set of constraints measures the difference 

between actual inventories of type (i,k) people and 

the requirement p.^ . The objective is selected so that the 

solution to the linear program will have e^ positive and 

dik equal to zero if w^y^ exceed p^ and e^ equal to 

zero and d^ positive if w^^y^ is less than p^ • The 

e . . and d., thus measure the excess and deficit in the 
lk lk 

type (i,k) people account. 

The variables m. . and n. . play a similar excess 
and deficit role in the type (i,j) billet account. 

The first set of constraints for the LP-P program is 



(A- 3) 



K 

x . . - T g, .w . , y, =0 for all ( i , j ) 

il kj lk^k ,J 



This again is a conservation constraint. There are 
w ik-^k *"YP e (i'k) people and gj\ of them are assigned to 
job j. Thus these constraints assure that the assignments 
to billets, x^ j , are consistent with actual manpower available. 

The second and third sets of constraints in LP-P 
are identical with the second and third sets of constraints 
in LP-B. They define the excess and deficit in people and 
billet accounts. 
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The form of the conservation constraint in LP-P allow 
for some simplifications. Notice that if all the y^. are 
nonnegative then (A-3) implies that all the x„ will be 
nonnegative. Thus we can use (A-3) to eliminate from 

the problem. With this simplification the constraints of 
LP-P are shown in Table A-2 . 

The objective for both linear programs will typically 
be a compromise between minimizing costs and minimizing the 
penalty that measures our departure from people and billet 
targets. The cost objective is 



(A- 4 ) 



K 

I 

k=l 



c k*k 



The penalty measure is 



(A- 5) 



where 



(A-6) 



l 

i=l j 


J 7T..n.. + y . .m. . 

13 13 13 13 


(i) 


ir.. = (0..b..) 1 

13 13 13 


(ii) 


y . . = (tj; . .b . . ) 
13 13 13 


(iii) 


Y ik - ^ik^ik^ 


(iv) 


°ik = ({ ikPik rl 



I I 
1 ' 

i=l k=l 



E l Y ik e ik + a ” d 



ik ik 
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Recall that 6, <p , and 5 are the percentage errors that 
lead to a unit loss. 

A balance between these two objectives is obtained by 
taking a weighted combination. 



K 

(A-7) A l c, y, + (1 - A) 
k=l K X 



K 



x ‘Ji + p ij m ij + Y ik e ik + °ik d ik } 



The problem LP-B has I x (J + 2K) constraints and 
K + I x (2K + 3J) variables. For our example with 1=6, 

J = 7 and K = 5 this means 102 constraints and 191 variables. 

The problem LP-P has I x (J + K) constraints and K + I x (2K + 2j) 
variables . For our example this works out to 72 constraints 
and 149 variables. LP-P is particularly easy to solve since 
one can always get a reasonable first basic solution by guessing 



the y^ and then choosing the d^ k , e^, m^j , and n . to 



il 



satisfy the constraints. 

The duals of both LP-B and LP-P appear to be easier to 
solve. The dual variables for each problem and each set of 
constraints is shown in Tables (A-l) and (A-2) respectively. 
The dual program for LP-B is 
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(A-8) maximize 



I K 



I J 



l l v ik p ik + l i 

i=i k=i 1K 1K i=i j=i 13 13 



K 



subject to 



£ u ik f kj 



for all (i,j) 



I I 

• i, u ik w ik + I, v ik u ik i Xc k for a11 k 

1=1 i=l 

- (1 - X) Y ik £ v ik £ (1 - A) a ik for all (i,k) 

- (l-A)ir^j £ q i j £ (1 - A) for all (i,j) 

Dual of LP-B 



The dual of LP-B has K + I x J constraints, I x (2K + J) 
variables, and I x (K + J) of these variables have upper and 
lower bounds. In our example this works out to 47 constraints, 
102 variables, and 72 variables subject to upper and lower bound 
constraints . 

The dual of LP-P is shown below. 
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(A-9 ) 



maximize 



I K 

y y v.,p., + 

i=i k=i lk lk 



subject to l v.-w.fc. + l 
i=l lK 1K i=l 



' (1 ' A) ^ik < V ik 



- (1 - X) TT. . < q . 

n — M i 



id 



Dual of LP-P 



I J 



l l q- b.. 

i=l j=l ^ x 3 

J 

^ij^kj w ik — ^ c k 



< (1 - X)a. k 



< (1 - X) y. . 
— ij 



for all k 

for all (k,i) 
for all (i,j) 



This problem has K constraints and I x (J + K) variables with 
upper and lower bounds. In our example, this would be 5 con- 
straints and 72 variables with upper and lower bounds. This 
problem may be quite easy to solve. 
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APPENDIX B 



UNCONSTRAINED QUADRATIC MODELS 



This appendix outlines two methods for selecting an allo- 
cation using the quadratic penalty functions. The two models are 
similar to the two linear programming models; they use the billet 
sharing or people sharing rules to go from an aggregate problem 
to a detailed allocation. 

The models do not have any inequality contraints. This 
is to insure that the solution can be obtained quickly. A full 
model with inequality constraints might appear to be more appro- 
priate, however, we must recall our ultimate objective is to cal- 
culate allocations in a relatively simple way. The model is not 
built on exact premises and it does not use precise data. It is there- 
fore not terribly important to be exact in choosing an allocation. 

We hope to have an easy and consistent way of choosing allocations 
so we can compare the effects of changing assumptions on the 
allocations . 

Both models calculate the billet assignments x^ and 
the accession rules y^. The actual allocation is determined by 
the billet sharing rule in model UQ-B 
Rule B. 



(B-l) 




f 



, .x. . . 

kj 13 
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In model UQ-P, the people sharing rule is used 
Rule P. 



(B-2 ) 



i 





These rules are discussed in Appendix A. 

The objective in these quadratic models is based on 
symmetric penalties; i.e. 0^ = and = ^ik* The 

objective is a combination of two terms: the cost term and 

a penalty term 



In the UQ-B model we minimize (B-3) subject to the conservation 
constants . 



K 



(B-3) 1 l c k y k + (1 - A) 




; j . .2 j K 





where 



(B-4) 




-2 




-2 




J 



(B-5) 
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The people sharing model, UQ-P, minimizes (B-3) subject 



to the constraints 



(B-6 ) 




K 

^j w ik^k ^ 



In UQ-P we can substitute (B-6) directly into (B-3) and thus 
solve a completely unconstrained problem. In UQ-B we minimize 
(B-3) subject to (B-5). 
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APPENDIX C 



SAMPLE SOLUTION 

This section contains some sample solutions using the 

unconstrained minimization described in Appendix B. 

The data, b^j, W ik' ^ik' are described in the paper. 

The p., and g?; , are calculated in the manner suggested in 
IK. K 

Section 12. 

The solutions are based on data that is largely subjective 
and does not, in any way, pretend to capture the situation that 
exists in the Navy. The intent of this section is to demonstrate 
the feasibility of the scheme. 

The target errors are expressed in percentage terms. The 
10,000% target error indicates that we do not care very much about 
meeting billit or personnel targets in the first two stages. 



Job 

Stage' Vs ^s Ss ^ 


1 


2 


3 


4 


5 


6 


7 


ENS 


10000 


10000 


10000 


10000 


10000 


10000 


10000 


LTJG 


10000 


10000 


10000 


10000 


10000 


10000 


10000 


LT 


10 


8 


6 


4 


4 


5 


6 


LTCDR 


8 


6.4 


4.8 


3.2 


3.2 


4 


4.8 


CDR 


7 


5.6 


4.2 


2.8 


2.8 


3.5 


4.2 


CAPT 


6 


4.8 


3.6 


2.4 


2.4 


3 


3.6 


Table C-l 


The 


inputs 


in % ; for billet requirements 



Job 

Stage 1 


2 


3 


4 


5 


ENS 0 


10000 


10000 


10000 


10000 


LTJG 10000 


10000 


10000 


10000 


10000 


LT 

LTCDR ^ 


7.5 


4.5 


5 


6 


CDR 8 


6 


3 . 6 


4 


4.8 


CAPT 7 


5.25 


3.15 


3.5 


4.2 


6 


4.5 


2.7 


3 


3 . 6 


Table C-2/ The 


inputs 


<}> ik in % ; 


for manpower 


inventorie 
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The billet assignments and the actual percentage errors 



using the people share rule are 



STAGE/ JOB 


1 


2 


3 


4 


5 


6 


7 


ENS 


1444 


0 


2231 


1000 


1791 


972 


0 


LTJG 


1498 


0 


1989 


895 


1669 


908 


1 


LT 


2164 


383 


2380 


19 76 


2021 


992 


399 


LTCDR 


1242 


444 


1246 


586 


977 


39 5 


452 


CDR 


1108 


482 


1000 


489 


578 


63 


884 


CAPT 


540 


529 


308 


150 


0 


0 


510 




Table C- 


■3. The billet assignment, 


x . 
13 






Table C-4. The percentage error in meeting billet requirements 



The accessions to the five classes are: 
Class 1 2345 

J 645 1162 5Z4 926 508 

Table C-5. Annual rate of accessions, 



These accessions produce the actual personnel inventories 
and the precentage error in meeting inventories are given below. 



Class 


1 


2 


3 


4 


5 


ENS 


1274 


2296 


1035 


1829 


1003 


LTJG 


1192 


2149 


969 


1712 


938 


LT 


1613 


2906 


1311 


2316 


1269 


LTCDR 


909 


1637 


738 


1304 


715 


CDR 


789 


1421 


641 


1132 


621 


CAPT 


349 


629 


284 


501 


275 



Table C-6. Personnel inventory by class and stage, 




Table C-7. Percentage errors in meeting personnel inventory budgets. 



The solution using the billet sharing rule produced similar 
results . The billet assignments and % error in billet assignments 
were 
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STAGE/JOB 

ENS 

LTJG 

LT 

LTCDR 

CDR 

CAPT 



1801 

1685 

2280 

1284 

1115 

493 



0 

0 

405 

370 

496 

554 



2074 

1941 

2423 

1294 

1036 

291 



940 

879 

1109 

596 

483 

147 



1388 

1298 

1399 

676 

101 

0 



1015 

949 

1021 

493 

64 

0 



0 

1 

499 

432 

1173 

492 



Table C-8 



STAGE /JOB 
ENS 
LTJG 
LT 

LTCDR 

CDR 

CAPT 




Table C-9 



The accessions using billet sharing were 

Class 12345 

638 1187 530 748 550 

Table C-10 
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These accession rates give rise to the personnel 



inventory and percentage error in personnel levels 



Class 


1 


2 


3 


4 


5 


ENS 


1261 


2345 


1048 


1478 


1087 


LTJG 


1180 


2194 


981 


1383 


1017 


LT 


1596 


2968 


1326 


1871 


1376 


LTCDR 


899 


1672 


747 


1054 


775 


CDR 


780 


1451 


649 


915 


673 


CAPT 


345 


642 


287 


405 


29 8 






Table 


C-ll 







Class 


1 


2 


3 


4 


5 


ENS 


“340 


9 


47 


51 


“l02 


LTJG 


67 


8 


23 


33 


19 


LT 


26 


“l 7 


29 


54 


23 


LTCDR 


14 


10 


1 


31 


24 


CDR 


17 


2 


2 


20 


“20 



Table C-12 



The analysis of this output should be directed toward 
constructive changes in the input data. Can we shift some 
billet assignments from LT to LTJG and LTCDR? Should we 
tighten up more on the pilot inventory? Should we change 
the sharing rules? Is it possible to alter the survivor 
fractions, and thereby improve the solution? 
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We see that there are a number of potential questions 
we can answer. The model is flexible and the calculations 
are rapid. It should be an excellent tool for analyzing 
manpower policy. 
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APPENDIX D 



ORGANIZATION OF PROGRAMS 



This appendix shows how the programs and files 
for the optimization problem are set up. 
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FILE STRUCTURE 



Component 

Number 


Dimension 


Symbol 


Description 


1 


(K,I, J) 


G 


People sharing array. +/G is 

an (I,K) element matrix of ones. 

. . i 

G is nonnegative: g^ is the 

fraction of type (i,k) people 

assigned to job(i,j) 


2 


(1, J) 


B 


j 

The billet requirements. b^_. 
is the number of type(i,j) jobs 
to be filled. 


3 


(1, J) 


0 


The range for a unit error in 

overfilling billets. If 

x. . = (1 + 0.01 . .)b. . there 
i] 13 13 

is a unit error. 0. . >0. 

13 

Expressed as %. 


4 


(I , K) 




The range for a unit error in 
overfilling personnel inventories 

z ik = (1 l°-°^ik )p ik a unit 
error is counted. <J>^ > 0. 

Expressed as % . 


5 


(K , I , J) 


F 


The job sharing array. F is 

nonnegative and +/[1]F is 

an (I,J) matrix of ones: f?; . 

3 

is the fraction of type ( i , j ) 
jobs to be filled by type(i,k) 
people . 


6 


( I » K) 


P 


The desired manpower inventory: 
p^ is the number type(i,k) 
people desired. 
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FILE STRUCTURE CONT. 



Component 

Number 


Dimens ional 


Symbol 


Description 


7 


(N,K) 


a 


The survivor fractions for 
each of the K classes . N 
is the maximum LOS. 


8 


I 


s 


The last year in each stage. 
Increasing positive numbers 
with S = N. 


9 


(I,K) 


W 


The element w., is the number 

lk 

of years a class k input expects 
to spend in stage i. 


10 


K 


y 


The calculated accessions rate; 
y^ is the number of accessions 
in manpower class k. 


11 


(I,K) 


Z 


The calculated inventory of 
each type of person. There are 
z ik type People. 


12 


(K,I,J) 


A 


The calculated allocation: a, 1 . 

kD 

is the number of type (i,k) 
people filling type(i,j) jobs. 


13 


(I,J) 


X 


The calculated billets filled. 
There are x^ people in job(i,j) . 


14 


(K , I , J) 


T 


The target allocation. We want 
t£j type (i,k) people in type(i,j) 
jobs . If things are consistent 
+/T equals P and +/[1]T equals B. 


15 


(I, J) 


P 


The range for a unit error in 

under filling billets. If 

x. . = (1-iii. . )b. . then an error of 1 
iD ID iD 

is counted . \b. . > 0 . 

ID “ 
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FILE STRUCTRE CONT . 



Component 

Number 


Dimensional 


Symbol 


Description 


16 


( I / K) 


6 


The range for a unit error in 
underfilling personnel inven- 
tories. I£ z. k = d-« lk )P lk 

then a unit error is counted. 
6 ik > °. 


17 


(K , I , J) 


C 


The element c, 1 . is the total 
cost of having a type (i,k) 
person in job(i,j) for one 
year . 
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FUNCTIONS 



Function 


Uses 


Computes 


Syntax and Description 


BILQUAD 


B, 0, <f>, F, 
P, W 


Y, X, Z, A 


BILQUAD ' FILENAMS 1 
Calculates "optimal" acces- 
sions y, personnel inventory 
(z), and allocation (A), 
given the data, objective, 
and the billet share rule (F) . 


BILSIM 


F,Y,X 


Z , A 


BILSIM' FILENAMS' 

Given Y and X (which pre- 
sumably satisfy the conser- 
vation constraint) , calculates 
allocation (A) and personnel 
inventory (Z), using the 
billet sharing rule. 


PGREC 


P/ G 


B, F, T 


PGREC 'FILENAMS' 
reconciles the file elements 
B, G, and T with P and G 


BFREC 


B, F 


P, G, T 


BFREC ' FILENAMS ' 
reconciles the file elements 
P, G, and T with the file 
elements B and F in the 
same file. 


PEOQUAD 


G/ B , 0 , (p r 

F, P, W 


Y, X, Z , A 


PEOQUAD 'FILENAMS' 
calculates optimal acces- 
sions (y) , personnel in- 
ventory (z), billet staffing 
(x) , and allocation (A) , 
given the data, and use of 
the people share rule . 
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FUNCTIONS CONT. 



FUNCTION 


USES 

COMPONENTS 


CALCULATES 

COMPONENT 


SYNTAX AND DESCRIPTION 


PEOSIM 


G, Y, W 


Z, X, A 


PEOSIM' FILENAMS ' 
calculates inventory z, 
staffing x and allocation 
A, where accessions and 
people share rule are given . 


WAITS 


a, S 


W 


WAITS ' FILENAMS ' 

Given the survivor frac- 
tions (a) and stage defini- 
tions (s), calculates w 
the expected waiting time 
in each state. 


PEOPRONT 


P, z 


lOOx(D-z) tP 


PEOPRONT ' FILENAMS ' 
calculates the percentage 
error in inventory. Compare 
with . 


BILPRONT 


B, X 


100x (B-x) vB 


BILPRONT 1 FILENAMS ' 
calculates the percentage 
error in meeting billet 
requirements. Compare 
with 9. 


ROUND 






ROUND XXX 

takes any array and rounds 
elements to integers. 
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